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1 Introduction 

Quantum chromodynamics, or QCD, as it is normally known in high energy 
physics, is the quantum held theory that describes the strong interactions. 
It is the SU(3 ) gauge theory of the current Standard Model for elementary 
particles and forces, SU(3) x SU(2)l x 17(1), which encompasses the strong, 
electromagnetic and weak interactions. The symmetry group of QCD, with 
its eight conserved charges, is referred to as color SU{3). As is characteristic 
of quantum held theories, each held may be described in terms of quantum 
waves or particles. 

Because it is a gauge held theory, the helds that carry the forces of QCD 
transform as vectors under the Lorentz group. Corresponding to these vector 
holds are the particles called gluons , which carry an intrinsic angular momen¬ 
tum, or spin, of one in units of h. The strong interactions are understood as 
the cumulative effects of gluons, interacting among themselves and with the 
quarks, the spin-1/2 particles of the Dirac quark helds. 

There are six quark helds of varying masses in QCD. Of these, three 
are called light quarks, in a sense to be defined below, and three heavy. 
The light quarks are the up (u), down (d) and strange (s), while the heavy 
quarks are the charm (c), bottom (6) and top (t). Their electric charge 
are famously e/ = 2e/3 (u,c,t) and e/ = —e/3 (d,s,b), with e the positron 

1 Contribution to appear in the Elsevier Encyclopedia of Mathematical Physics. 
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charge. The gluons interact with each quark held in an identical fashion, 
and the relatively light masses of three of the quarks provide the theory with 
a number of approximate global symmetries that profoundly influence the 
manner in which QCD manifests itself in the Standard Model. 

These quark and gluon fields and their corresponding particles are enu¬ 
merated with complete confidence by the community of high energy physi¬ 
cists. Yet none of these particles has ever been observed in isolation, as 
one might observe a photon or an electron. Rather, all known strongly in¬ 
teracting particles are colorless; most are mesons, combinations with the 
quantum numbers of a quark q and a antiquark q', or baryons with the quan¬ 
tum numbers of (possibly distinct) combinations of three quarks qq'q". This 
feature of QCD, that its underlying fields never appear as asymptotic states, 
is called confinement. The very existence of confinement required new ways 
of thinking about held theory, and only with these was the discovery and 
development of QCD possible. 

2 The Background of QCD 

The strong interactions have been recognized as a separate force of nature 
since the discovery of the neutron as a constituent of atomic nuclei, along 
with the proton. Neutrons and protons (collectively, nucleons) possess a 
force, attractive at intermediate distances and so strong that it overcomes the 
electric repulsion of the protons, each with charge e. A sense of the relative 
strengths of the electromagnetic and strong interactions may be inferred from 
the typical distance between mutually-repulsive electrons in an atom, about 
10~ 8 cm, and the typical distance between protons in a nucleus, of order 
10“ 13 cm. 

The history that led up to the discovery of QCD is a fascinating one, 
beginning with Yukawa’s 1935 theory of pion exchange as the source of the 
forces that bind nuclei, still a useful tool for low energy scattering. Other 
turning points include the creation of nonabelian gauge theories by Yang and 
Mills in 1954, the discovery of the quantum number known as strangeness, 
the consequent development of the quark model and then the proposal of 
color as a global symmetry. The role of point-like constituents in hadrons 
was foreshadowed by the identification of electromagnetic and weak currents 
and the analysis of their quantum-mechanical algebras. Finally, the obser¬ 
vation of scaling in deep-inelastic scattering, which we will describe below, 
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made quantum chromodynamics, with color as a local symmetry, the unique 
explanation of the strong interactions, through its property of asymptotic 
freedom. 


3 The Lagrangian and its Symmetries 

The QCD Lagrangian may be written as 

n f i \ 

C =12qf ( iJ P[ A \ -™/) 9 /-xTr[F^(hL)]-- (. B a (A)f+c b 

f= i z z 

( 1 ) 

with lp[A} = 7 • d + igg'y ■ A the covariant derivative in QCD. The are the 
Dirac matrices, satisfying the anticonnnutation relations, [ 7 ^, 7 ^]+ = 2g fiu . 
The 577(3) gluon fields are A M 7L^T a , where T a are the generators 

of SU( 3) in the fundamental representation. The held strengths F^ v [A] = 
d^A v — d v A M + ig^A^, A u \ specify the three- and four-point gluon couplings 
of nonabclian gauge theory. In QCD, there are rif = 6 flavors of quark fields, 
qj, with conjugate qf = qj'y 0 . 

The first two terms in the expression m make up the classical Lagrangian, 
followed by the gauge fixing term, specified by a (usually, but not necessarily 
linear) function B a (A ), and the ghost Lagrangian. The ghost (antighost) 
fields c a (c a ) carry the same adjoint index as the gauge fields. 

The classical QCD Lagrangian before gauge fixing is invariant under the 
local gauge transformations 



= —du^i(x) D x (x) 

+ D(x) A!^{x) Ft 1 (x) 
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= A^ix) - d^8oi(x) 

+ ig s [Sa^x)^^)] + .. 

$0*0 

= VL{x)ij ipj{x) = -0 

i(x) + ig s 8a(x)ij ^j{x) + 

Sa(x ) 

= $ a a(x) T a . 

a= 1 



The full QCD action including gauge fixing and ghost terms is also invariant 
under the BRST transformations, with an anticommuting variable, 

8A^ a (fiabd/j, T gA^ c fabc) Q) 

8c a 2 9Cube c c S £ Sc A.B a 

Hi = W [Tblij cHj , (3) 
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and with f abc the SU( 3) structure constants. The Jacobean of these trans¬ 
formation is unity. 

In addition, neglecting masses of the light quarks, u, d and s, the QCD 
Lagranian has a class of global flavor and chiral symmetries, the latter 
connecting left- and right-handed components of the quark fields, i(> l,r = 

(1/2) (1 T 75 ) V’, 

= e ia ^ x ), P = 0,1. (4) 

Here power P = 0 describes phase, and P = 1 chiral, transformations. 
Both transformations can be extended to transformations among the light 
flavors, by letting ^ become a vector, and a an element in the Lie Algebra 
of SU(M), with M — 2 if we take only the u and d quarks, and M = 3 if we 
include the somewhat heavier strange quark. These symmetries, not to be 
confused with the local symmetries of the Standard Model, are strong isospin 
and its extension to the “eightfold way”, which evolved into the (3-)quark 
model of Gell-Mann and Zweig. The many successes of these formalisms are 
automatically incorporated into QCD. 

4 Green Functions, Phases and Gauge Invari¬ 
ance 

In large part, the business of quantum held theory is to calculate Green 
functions, 

G n (xi... x n ) = (0 | T ($i(xi)... ... & n (x n ) ) | 0) , (5) 

where T denotes time ordering. The $j(x) are elementary fields, such as 
A or qf , or composite holds, such as currents like Such a 

Green function generates amplitudes for the scattering of particles of def ini te 
momenta and spin, when in the limit of large times the ^-dependence of the 
Green function is that of a plane wave. For example, we may have in the 
limit x® —> 00 , 

G n (sq ... x n ) * 

0i(p, A) e ipXi ((p, A) | T ($i(aq)... <I>i_i(aq_i) $ i+ i(ic i+ i)... $ n (x n ) ) | 0) , 

( 6 ) 

where (f>i(p, A) is a solution to the free-held equation for held <Fj, character¬ 
ized by momentum p and spin A. (An inegral over possible momenta p is 
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understood.) When this happens for field i, the vacuum state is replaced 
by \(p, A)), a particle state with precisely this momentum and spin; when it 
occurs for all fields, we derive a scattering (5) matrix amplitude. In essence, 
the statement of confinement is that Green functions with fields qj(x) never 
behave as plane waves at large times in the past or future. Only Green func¬ 
tions of color singlet composite fields, invariant under gauge transformations, 
are associated with plane wave behavior at large times. 

Green functions remain invariant under the BRST transformations ©, 
and this invariance implies a set of Ward identities 

6 n 

£ (0 | T (^(Xi) . . . S B RS*i(Zi) • • • $n(Xn) ) | 0) = 0 . (7) 

<W) ^ 


The variation of the antighost as in © is equivalent to an infinitesimal change 
in the gauge fixing term; variations in the remaining fields all cancel single¬ 
particle plane wave behavior in the corresponding Green functions. These 
identities then ensure the gauge-invariance of the perturbative S-matrix, a 
result that turns out to be useful despite confinement. 

To go beyond a purely perturbative description of QCD, it is useful to 
introduce a set of nonlocal operators that are variously called nonabelian 
phases, ordered exponentials and Wilson lines, 


U c (z,y) = P ex p 


-ig s 



dx M A^(x) 


J 


( 8 ) 


where C is some self-avoiding curve between y and 0 . The LPs transform at 
each end linearly in nonabelian gauge transformations fl(x) at that point, 


U' c (z,y) = Sl(z)Uc(z,y)tt '(y). 


( 9 ) 


Especially interesting are closed curves C, for which z = y. The phases 
about such closed loops are, like their abelian counterparts, sensitive to the 
magnetic flux that they enclose, even when the field strengths vanish on the 
curve. 


5 QCD at the Shortest and Longest Distances 

Much of the fascination of QCD is its extraordinary variation of behavior 
at differing distance scales. Its discovery is linked to asymptotic freedom, 
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which characterizes the theory at the shortest scales. Asymptotic freedom 
also suggests (and in part provides) a bridge to longer distances. 

Most analyses in QCD begin with a path integral formulation in terms of 
the elementary fields <F a = qy ..., 


G n (x i ,(z j ,y j )) = j 


n 


a=q,q,G,c,c 


n*ife)nc Cj fe,®)e iSa “. 

i j 


( 10 ) 


with S'qcd the action. Perturbation theory keeps only the kinetic Lagrangian, 
quadratic in fields, in the exponent, and expands the potential terms in 
the coupling. This procedure produces Feynman diagrams, with vertices 
corresponding to the cubic and quartic terms in the QCD Lagrangian Q. 

Most nonperturbative analyses of QCD require studying the theory on a 
Eucliean, rather than Minkowski space, related by an analytic continuation 
in the times x°,y°,z° in G n from real to imaginary values. In Euclidean 
space, we find, for example, classical solutions to the equations of motion, 
known as instantons, that provide nonperturbative contributions to the path 
integral. Perhaps the most flexible nonperturbative approach approximates 
the action and the measure at a lattice of points in four dimensional space. 
For this purpose, integrals over the gauge fields are replaced by averages over 
“gauge links”, of the form of Eq. (JHJ) between neighboring points. 

Perturbation theory is most useful for processes that occur over short 
time scales and at high relative energies. Lattice QCD, on the other hand, 
can simulate processes that take much longer times, but is less useful when 
large momentum transfers are involved. The gap between the two methods 
remains quite wide, but between the two they have covered enormous ground, 
enough to more than confirm QCD as the theory of strong interactions. 


5.1 Asymptotic freedom 

Quantum chromodynamics is a renormalizable field theory, which implies 
that the coupling constant g must be defined by its value at a renormalization 
scale , and is denoted g(/i). Usually, the magnitude of ct s (/i) = g 2 /4:7T, is 
quoted at p = mz, where it is approximately 0.12. In effect, g(y) controls 
the amplitude that connects any state to another state with one more or one 
fewer gluon, including quantum corrections that occur over time scales from 
zero up toh/y (if we measure /i in units of energy). The QCD Ward identities 
mentioned above ensure that the coupling is the same for both quarks and 
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gluons, and indeed remains the same in all terms in the Lagrangian, ensuring 
that the symmetries of QCD are not destroyed by renormalization. 

Quantum corrections to gluon emission are not generally computable di¬ 
rectly in renormalizable theories, but their dependence on /i is computable, 
and is a power series in a s (n) itself, 


A 


> doi s (/r) 
dji 2 


= -h 




4:71 


-h 


(47t)- 


+ ... = /3(a s ), 


( 11 ) 


where bo = 11 — 2n//3 and b\ = 2(31 — 19nj/3). The celebrated minus signs 
on the right-hand side are associated with both the spin and self-interact ions 
of the gluons. 

The solution to this equation provides an expression for a s at any scale 
(ii in terms of its value at any other scale fi 0 . Keeping only the lowest-order, 
b 0 , term, we have 


l) 


(ho) 

1 + (& 0 /4vr) In (j4/f4) 


4n 

b 0 ln(/jy/AQ CD ) 


( 12 ) 


where in the second form, we have introduced Aqcd, the scale parameter 
of the theory, which embodies the condition that we get the same coupling 
at scale Hi no matter which scale /io we start from. Asymptotic freedom 
consists of the observation that at larger renormalization masses /i, or cor¬ 
respondingly shorter time scales, the coupling weakens, and indeed vanishes 
in the limit /i —> oo. The other side of the coin is that over longer times or 
lower momenta the coupling grows. Eventually, near the pole at = Aqcd, 
the lowest-order approximation to the running fails, and the theory becomes 
essentially nonperturbative. Thus the discovery of asymptotic freedom sug¬ 
gested, although it certainly doesn’t prove, that QCD is capable of producing 
very strong forces, and confinement at long distances. Current estimates of 
Aq CD are around 200 MeV. 


5.2 Spontaneous breaking of chiral symmetry 

The number of quarks and their masses is an external input to QCD. In the 
Standard Model masses are provided by the Higgs mechanism, but in QCD 
they are simply parameters. Because the Standard Model has chosen several 
of the quarks to be especially light, QCD incorporates the chiral symmetries 
implied by Eq. Q (with P = 1). In the limit of zero quark masses, these 
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symmetries becomes exact, respected to all orders of perturbation theory, 
that is, for any finite number of gluons emitted or absorbed. 

At distances on the order to 1 /Aqcd, however, QCD cannot respect chiral 
symmetry, which would require each state to have a degenerate partner with 
the opposite parity, something not seen in nature. Rather, QCD produces, 
nonperturbatively, nonzero values for matrix elements that mix right- and 
left-handed fields, such as (0 |wlWr| 0), with u the up-quark held. Pions are 
the Goldstone bosons of this symmetry, and may be thought of as ripples in 
the chiral condensate, rotating it locally as they pass along. The observation 
that these Goldstone bosons are not exactly massless is due to the current 
masses of the quarks, their values in £qcd- The (chiral perturbation the¬ 
ory) expansion in these light quark masses also enables us to estimate them 
quantitatively: 1.5 < m u < 4 MeV, 4 < < 8 MeV, 80 < m s < 155 

MeV. These are the light quarks, with masses smaller than Aqcd- (Like a s , 
the masses are renormalized; these are quoted from Eidelman (2004) with 
/i = 2 GeV.) For comparison, the heavy quarks have masses m c ~ 1 — 1.5 
GeV, nib ~ 4 — 4.5 GeV, and rn t ~ 180 GeV (the giant among the known 
elementary particles). 

Although the mechanism of the chiral condensate (and in general other 
nonperturbative aspects of QCD) has not yet been demonstrated from first 
principles, a very satisfactory description of the origin of the condensate, and 
indeed of much hadronic structure has been given in terms of the attractive 
forces between quarks provided by instantons. The actions of instanton so¬ 
lutions provide a dependence exp[—87r 2 /g 2 ] in Euclidean path integrals, and 
so are characteristically nonperturbative. 

5.3 Mechanisms of confinement 

As described above, confinement is the absence of asymptotic states that 
transform nontrivially under color. The full spectrum of QCD, however, is a 
complex thing to study, and so the problem has been approached somewhat 
indirectly. A difficulty is the same light quark masses associated with approx¬ 
imate chiral symmetry. Because the masses of the light quarks are far below 
the scale Aqcd at which the perturbative coupling blows up, light quarks are 
created freely from the vacuum and the process of “hadronization”, by which 
quarks and gluons form mesons and baryons, is both nonperturbative and 
relativistic. It is therefore difficult to approach in both perturbation theory 
and lattice simulations. 



Tests and studies of confinement are thus normally formulated in trun¬ 
cations of QCD, typically with no light quarks. The question is then refor¬ 
mulated in a way that is somewhat more tractable, without relativistic light 
quarks popping in and out of the vacuum all the time. In the limit that 
its mass becomes infinite compared to the natural scale of fluctuations in 
the QCD vacuum, the propagator of a quark becomes identical to a phase 
operator, (JSJ) , with a path C corresponding to a constant velocity. This ob¬ 
servation suggests a number of tests for confinement that can be implemented 
in the lattice theory. The most intuitive is the vacuum expectation value of 
a “Wilson loop”, consisting of a rectangular path, with sides along the time 
direction, corresponding to a heavy quark and antiquark at rest a distance R 
apart, and closed at some starting and ending times with straight lines. The 
vacuum expectation value of the loop then turns out to be the exponential of 
the potential energy between the quark pair, multiplied by the elapsed time, 



(13) 


When V(R) oc R (“area law” behavior), there is a linearly-rising, confining 
potential. This behavior, not yet proven analytically yet well-confirmed on 
the lattice, has an appealing interpretation as the energy of a “string”, con¬ 
necting the quark and antiquark, whose energy is proportional to its length. 

Motivation for such a string picture was also found from the hadron spec¬ 
trum itself, before of any of the heavy quarks were known, and even before 
the discovery of QCD, from the observation that many mesonic {qq') states 
lie along “Regge trajectories”, which consist of sets of states of spin J and 
mass 77ij that obey a relation 


J = a'm 2 j , 


(14) 


for some constant a '. Such a relation can be modeled by two light particles 
(“quarks”) revolving around each other at some constant (for simplicity fixed 
nonrelativistic) velocity v 0 and distance 2 R, connected by a “string” whose 
energy per unit length is a constant p. 

Suppose the center of the string is stationary, so the overall system is at 
rest. Then neglecting the masses, the total energy of the system is M — 2 Rp. 
Meanwhile, the momentum density per unit length at distance r from the 
center is v(r) = (r/R)vo, and the total angular momentum of the system is 



( 15 ) 
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and for such a system, m is indeed satisfied. Quantized values of angular 
momentum J give quantized masses mj, and we might take this as a sort 
of “Bohr model” for a meson. Indeed, string theory has its origin in related 
consideration in the strong interactions. 

Lattice data are unequivocal on the linearly rising potential, but it re¬ 
quires further analysis to take a lattice result and determine what field con¬ 
figurations, string-like or not, gave that result. Probably the most widely 
accepted explanation is in terms of an analogy to the Meissner effect in su¬ 
perconductivity, in which type 11 superconductors isolate magnetic flux in 
quantized tubes, the result of the formation of a condensate of Cooper pairs 
of electrons. If the strings of QCD are to be made of the gauge field, they must 
be electric (F M °) in nature to couple to quarks, so the analogy postulates a 
“dual” Meissner effect, in which electric flux is isolated as the result of a con¬ 
densate of objects with magnetic charge (producing nonzero F y ). Although 
no proof of this mechanism has been provided yet, the role of magnetic fluc¬ 
tuations in confinement has been widely investigated in lattice simulations, 
with encouraging results. Of special interest are magnetic field configura¬ 
tions, monopoles or vortices, in the Z 3 center of SU( 3), exp[z7rfc/3] / 3X 3 , 
k — 0,1,2. Such configurations, even when localized, influence closed gauge 
loops m through the nonabelian Aharonov-Bohm effect. Eventually, of 
course, the role of light quarks must be crucial for any complete description 
of confinement in the real world, as emphasized by Gribov. 

Another related choice of closed loop is the “Polyakov loop”, implemented 
at finite temperature, for which the path integral is taken over periodic field 
configurations with period 1/T, where T is the temperature. In this case, 
the curve C extends from times t — 0 to t — 1/T at a fixed point in space. In 
this formulation it is possible to observe a phase transition from a confined 
phase, where the expectation is zero, to a deconhned phase, where it is 
nonzero. This phase transition is currently under intense experimental study 
in nuclear collisions. 


6 Using Asymptotic Freedom: Perturbative 
QCD 

It is not entirely obvious how to use asymptotic freedom in a theory that 
should (must) have confinement. Such applications of asymptotic freedom 
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go by the term perturbative QCD, which has many applications, not the 
least as a window to extensions of the Standard Model. 

6.1 Lepton annihilation and infrared safety 

The electromagnetic current, J fl — e fQfhj,Qf is a gauge invariant operator, 
and its correlation functions are not limited by confinement. Perhaps the 
simplest application of asymptotic freedom, yet of great physical relevance, 
is the scalar two-point function, 

<Q) = Y jd A xe-^- x (0\T(J^0)J fl (x))\0) . (16) 

The imaginary part of this function is related to the total cross section for 
the annihilation process e + e _ —» hadrons in the approximation that only 
one photon takes part in the reaction. The specific relation is ctqcd = 
(e 4 /Q 2 ) Ini7r(Q 2 ), which follows from the optical theorem, illustrated in Fig. 
□ The perturbative expansion of the function 7 r(Q) depends, in general, on 
the mass scales Q and the quark masses m/ as well as on the strong coupling 
a s (/x) and on the renormalization scale //. We may also worry about the 
influence of other, truly nonperturbative scales, proportional to powers of 
Aq CD . At large values of Q 2 , however, the situation simplifies greatly, and 
dependence on all scales below Q is suppressed by powers of Q. This may 
be expressed in terms of the operator product expansion, 

(0 | T (J»( 0) J^x)) | 0) = J2(x 2 )~ 3+dl/2 Cj(x V , Oi s (fi)) (0 | 0,(0) | 0) , (17) 

o, 

where dj is the mass dimension of operator O,, and where the dimension¬ 
less coefficient functions Cj incorporate quantum corrections. The sum over 
operators begins with the identity (dj = 0), whose coefficient function is 
identified with the sum of quantum corrections in the approximation of zero 
masses. The sum continues with quark mass corrections, which are sup¬ 
pressed by powers of at least m 2 /Q 2 , for those flavors with masses below Q. 
Any QCD quantity that has this property, remaining finite in perturbation 
theory when all particle masses are set to zero, is said to be infrared safe. 

The effects of quarks whose masses are above Q are included indirectly, 
through the couplings and masses observed at the lower scales. In summary, 
the leading power behavior of n(Q), and hence of the cross section, is a 
function of Q, fi and a s (/x) only. Higher order operators whose vacuum matrix 
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elements receive nonperturbative corrections include the gluon condensate, 
identified as the product a s (/d)G a pG al3 oc Aq CD . 

Once we have concluded that Q is the only physical scale in 7r, we may 
expect that the right choice of the renormalization scale is fi — Q. Any 
observable quantity is independent of the choice of renormalization scale, p, 
and neglecting quark masses, the chain rule gives 


da(Q/p, aq(/i)) _ da da 

V -—h 2 p(a s ) —— 

did dp oa s 


= 0 , 


(18) 


which shows that we can determine the beta function directly from the per¬ 
turbative expansion of the cross section. Defining a = a s (p)/ir, such a 
perturbative calculation gives 


Im 7 t(Q 2 ) 


3 

47T 



^1 ~\~ CL ~\h 0? 


^1.986 - 0.11571/ - 


(£» 0 /47 t) In 



,( 19 ) 


with b 0 as above. Now, choosing p — Q, we see that asymptotic freedom 
implies that when Q is large, the total cross section is given by the lowest 
order, plus small and calculable QCD corrections, a result that is borne 
out in experiment. Comparing experiment to an expression like (USD , one 
can measure the value of a s (Q), and hence, with Eq. (jT2j) . a s (p) for any 
p 3> Aqcd- Fig. [21 shows a recent compilation of values of a s from this kind 
of analysis in different experiments at different scales, clearly demonstrating 
asymptotic freedom. 


6.2 Factorization, scaling and parton distributions 

One step beyond vacuum matrix elements of currents are their expectation 
values in single particle states, and here we make contact with the discovery 
of QCD, through scaling. Such expectations are relevant to the class of ex¬ 
periments known as deep-inelastic scattering, in which a high energy electron 
exchanges a photon with a nucleon target. All QCD information is contained 
in the tensor matrix element 

WT(P,0) = ^ Y, / d 4 xe~ iq ' x (p,a\ J^O) J v {x) \p,a) , (20) 

with q the momentum transfer carried by the photon, and p, a the momentum 
and spin of the target nucleon, N. This matrix element is not infrared safe, 
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since in depends in principle on entire history of the nucleon state. Thus, it 
is not accessible to direct perturbative calculation. 

Nevertheless, when the scattering involves a large momentum transfer 
compared to Aqcd we may expect a quantum-mechanical incoherence be¬ 
tween the scattering reaction, which occurs (by the uncertainty principle) at 
short distances, and the forces that stabilize the nucleon. After all, we have 
seen that the latter, strong forces, should be associated with long distances. 
Such a separation of dynamics, called factorization, can be implemented in 
perturbation theory, and is assumed to be a property of full QCD. Fac¬ 
torization is illustrated schematically in Fig. |21 Of course, short- and long 
distances are relative concepts, and the separation requires the introduction 
of a so-called factorization scale, p F , not dissimilar to the renormalization 
scale described above. For many purposes, it is convenient to choose the two 
equal, although this is not required. 

The expression of factorization for deep-inelastic scattering is 

W& u (p,q)= J2 [ d £ c r(t,P,q,VF,as{p F )) fi /N (Z,V F ) , ( 21 ) 

where the functions C pi/ (the coefficient functions) can be computed as an 
expansion in a s (p F ), and describe the scattering of the partons , quarks and 
gluons, of which the target is made. The variable £ ranges from unity down 
to x = -q 2 /2p-q > 0, and has the interpretation of the fractional momentum 
of the proton carried by parton i. (Here — q 2 = Q 2 is positive.) The parton 
distributions fi /n can be defined in terms of matrix elements in the nucleon, 
in which the currents are replaced by quark (or antiquark or gluon) fields, as 

1 • -I- 

fq/N{x, n) = — J dx e~ zXxp (p, a \ q(Xn)U n (nX , 0) n • 7 q( 0) \p,a) . (22) 

n p is a light-like vector, and U n a phase operator whose path C is in the 
n-direction. The dependence of the parton distribution on the factorization 
scale is through the renormalization of the composite operator consisting of 
the quark fields, separated along the light-cone, and the nonabelian phase 
operator U n (n\, 0), which renders the matrix element gauge invariant by Eq. 

By combining the calculations of the C’s and data for W^'\ we can infer 
the parton distributions, fi/N- Important factorizations of a similar sort also 
apply to some exclusive processes, including amplitudes for elastic pion or 
nucelon scattering at large momentum transfer. 
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Eq. (1211) has a number of extraordinary consequences. First, because the 
coefficient function is an expansion in a s , it is natural to choose /i 2 F 
p-q (when x is of order unity). When Q is large, we may approximate C 1 -' 1 ' by 
its lowest order, which is first order in the electromagnetic coupling of quarks 
to photons, and zeroth order in a s . In this approximation, dependence on Q 
is entirely in the parton distributions. But such dependence is of necessity 
weak (again for x not so small as to produce another scale), because the 
fiF dependence of /i/jv(£, Hf) must be compensated by the pp dependence 
of Cf*', which is order a s . This means that the overall Q dependence of 
the tensor is weak for Q large when x is moderate. This is the scaling 
phenomenon that played such an important role in the discovery of QCD. 

6.3 Evolution: beyond scaling 

Another consequence of the factorization (EH), or equivalently of the oper¬ 
ator definition, EH is that the /i F -dependence of the coefficient functions 
and the parton distributions are linked. As in the lepton annihilation cross 
section, this may be thought of as due to the independence of the physically- 
observable tensor Wpj' from the choice of factorization and renormalization 
scales. This implies that the p,p- dependence of fi/N may be calculated per- 
turbatively since it must cancel the corresponding dependence in C % . The 
resulting relation is coventionally expressed in terms of the evolution equa¬ 
tions, 



(23) 


where P a b(0 are calculable as power series, now known up to a^. This rela¬ 
tion expands the applicability of QCD from scales where parton distributions 
can be inferred directly from experiment, to arbitrarily high scales, reach¬ 
able in accelerators under construction or in the imagination, or even on the 
cosmic level. 

At very high energy, however, the effective values of the variable x can 
become very small and introduce new scales, so that eventually the evolution 
of Eq. m fails- The study of nuclear collisions may provide a new high- 
density regime for quantum chromodynamics, which blurs the distinction 
between perturbative and nonperturbative dynamics. 
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6.4 Inclusive production 

Once we have evolution at our disposal, we can take yet another step, and 
replace electroweak currents with any operator from any extension of QCD, 
in the standard model or beyond, that couples quarks and gluons to the 
particles of as-yet unseen fields. Factorization can be extended to these 
situations as well, providing predictions for the production of new particles, 
F of mass M, in the form of factorized inclusive cross sections, 

CTAB^F(M)(M,Pa,Pb) = X f d£ad£bfi/A(€a,V) fj/B(£b, AO 

i,j=QfQf,G 

XH i:j ^ F (M)(XaPA, X b p B , M, /X, « s (/x)) , (24) 

where the functions H^p may be calculated perturbatively, while the f t /A 
and fj/B parton distributions are known from a combination of lower-energy 
observation and evolution. In this context, they are said to be universal , 
in that they are the same functions in hadron-hadron collisions as in the 
electron-hadron collisions of deep-inelastic scattering. In general the calcula¬ 
tion of hard-scattering functions H l3 is quite nontrivial beyond lowest order 
in a s . The exploration of methods to compute higher orders, currently as far 
as a,, has required extraordinary insight into the properties of multidimen¬ 
sional integrals. 

The factorization method helped predict the observation of the W and 
Z bosons of electroweak theory, and the discovery of the top quark. The 
extension of factorization from deep-inelastic scattering to hadron production 
is nontrivial; and indeed it only holds in the limit that the velocities, $, of 
the colliding particles approach the speed of light in the center-of-momentum 
frame of the produced particle. Corrections to the relation (ED are then at 
the level of powers of Pi — 1, which translates into inverse powers of the 
invariant mass(es) of the produced particle(s) M. Factorizations of this sort 
do not apply to low velocity collisions. Arguments for this result rely on 
relativistic causality and the uncertainty principle. The creation of the new 
state happens over time scales of order 1/M. Before that well-defined event, 
the colliding particles are approaching at nearly the speed of light, and hence 
cannot affect the distributions of each others’ partons. After the new particle 
is created, the fragments of the hadrons recede from each other, and the 
subsequent time development, when summed over all possible final states 
that include the heavy particle, is finite in perturbation theory as a direct 
result of the unitarity of QCD. 
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6.5 Structure of hadronic final states 


A wide range of semi-inclusive cross sections are defined by measuring prop¬ 
erties of final states that depend only on the flow of energy, and which bring 
QCD perturbation theory to the threshold of nonperturbative dynamics. 
Schematically, for a state N = \ki.-.k^), we define S(N ) = 
where s(f2) is some smooth function of directions. We generalize the e + e - 
annihilation case above, and define a cross section in terms of a related, but 
highly nonlocal, matrix element, 


da(Q) 

ds = a ° 


d 4 xe 


—iQ'X 


j m (o)<W / d 2 n s (n)S(n) -s) j^x) 


(25) 


where a 0 is a zeroth order cross sections, and where £ is an operator at 
spatial infinity, which measures the energy flow of any state in direction 12: 
£{y£)\ki ... /cat) = (1 /Q) /Cj°5 2 (12—12j). This may seem a little complicated, 

but like the total annihilation cross section, the only dimensional scale on 
which it depends is Q. The operator £ can be defined in a gauge invariant 
manner, through the energy-momentum tensor for example, and has a mean¬ 
ing independent of partonic final states. At the same time, this sort of cross 
section may be implemented easily in perturbation theory, and like the total 
annihilation cross section, it is infrared safe. To see why, notice that when 
a massless ( k 2 = 0) particle decays into two particles of momenta xk and 
(1 — x)k (0 < x < 1), the quantity S is unchanged, since the sum of the new 
energies is the same as the old. This makes the observable S(N) insensitive 
to processes at low momentum transfer. 

For the case of leptonic annihilation, the lowest-order perturbative contri¬ 
bution to energy flow requires no powers of a s , and consists of an oppositely- 
moving quark and antiquark pair. Any measure of energy flow that includes 
these configurations will dominate over correlations that require a s correc¬ 
tions. As a result, QCD predicts that in most leptonic annihilation events 
energy will flow in two back-to-back collimated sets of particles, known as 
jets. In this way, quarks and gluons are observed clearly, albeit indirectly. 

With varying choices of S, many properties of jets, such as their distri¬ 
butions in invariant mass, and the probabilities and angular distributions of 
multijet events, and even the energy dependence of their particle multiplici¬ 
ties, can be computed in QCD. This is in part because hadronization is dom¬ 
inated by the production of light quarks, whose production from the vacuum 
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requires very little momentum transfer. Paradoxically, the very lightness of 
the quarks is a boon to the use of perturbative methods. All these considera¬ 
tions can be extended to hadronic scattering, and jet and other semi-inclusive 
properties of final states also computed and compared to experiment. 

7 Conclusions 

Quantum chromodynamics is an extremely broad held, and this article has 
hardly scratched the surface. The relation of QCD-like theories to supersym¬ 
metric and string theories, and implications of the latter for confinement and 
the computation of higher order perturbative amplitudes, have been some 
of the most exciting developments of recent years. As another example, we 
note that the reduction of the heavy quark propagator to a nonabclian phase, 
noted in our discussion of confinement, is related to an additional symmetries 
of heavy quarks in QCD, with many consequences for the analysis of their 
bound states. Of the bibliography given below, one may mention the four 
volumes of Shifman (2001,2002), which communicate in one place a sense of 
the sweep of work in QCD. 

Our confidence in QCD as the correct description of the strong interac¬ 
tions is based on a wide variety of experimental and observational results. At 
each stage in the discovery, confirmation and exploration of QCD, the math¬ 
ematical analysis of relativistic quantum held theory entered new territory. 
As is the case for gravity or electromagnetism, this period of exploration is 
far from complete, and perhaps never will be. 
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Figure Captions 

Fig. 1 First line: schematic relation of lowest order e + e _ annihilation 
to sum over quarks q, each with electric charge e q . Second line pertur¬ 
bative unitarity for the current correlation function 7 r(Q). 

Fig. 2 Experimental variation of the strong coupling with scales; from 
Bethke (2004). 

Fig. 3 Schematic depiction of factorization in deep-inelastic scattering. 


See also 

AdS/CFT Correspondence. Aharanov-Bohm effect. BRST Quan¬ 
tization. Chiral symmetry. Current algebras. Dirac field and Dirac 
operator. Euclidean field theory. Gauge theory. Instantons in 
gauge theory. Lattice gauge theory. Operator product expansion. 
Perturbation theory and techniques. Perturbative renormalization 
and BRST. Renormalization, general theory. Scattering: funda¬ 
mental concepts and tools. Scattering: asymptotic completeness 
and bound states. Quantum field theory overview. Seiberg-Witten 
theory. Standard model of particle physics. 


19 



Key Words 

Quantum chromodynamics 

Gauge theory 

Asymptotic freedom 

Confinement 

Infrared safety 

Factorization 

Evolution 

Jets 

Renormalization 
Chiral symmetry 
Strong coupling 
High energy scattering 


20 



a(Q) = 

n(Q)= 5= 

m 


L 

q 


3 + 


hAAA; 


^ n(Q) 


C^m| = Im -O’" 

= lm( Q 

+ -< 2 >- + ■ ■ ■) 


Figure 1: 


21 









Figure 2: 


22 





















23 














